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Abstract 

Characterizing the global maximum of weighted sum-rate (WSR) for the i^-user Gaussian interference channel 
(GIC), with the interference treated as Gaussian noise, is a key problem in wireless communication. However, due 
to the users' mutual interference, this problem is in general non-convex and thus cannot be solved directly by 
conventional convex optimization techniques. In this paper, by jointly utilizing the monotonic optimization and rate 
profile techniques, we develop a new framework to obtain the globally optimal power control and/or beamforming 
solutions to the WSR maximization problems for the GICs with single-antenna transmitters and single-antenna 
receivers (SISO), single-antenna transmitters and multi-antenna receivers (SIMO), or multi-antenna transmitters 
and single-antenna receivers (MISO). It is assumed that the transmitted signals have circularly symmetric complex 
Gaussian distributions and are independent over time. Different from prior work, this paper proposes to maximize 
the WSR in the achievable rate region of the GIC directly by exploiting the facts that the achievable rate region 
is a "normal" set and the users' WSR is a "strictly increasing" function over the rate region. Consequently, the 
WSR maximization is shown to be in the form of monotonic optimization over a normal set and thus can be solved 
globally optimally by the existing outer polyblock approximation algorithm. However, an essential step in the 
algorithm hinges on how to efficiently characterize the intersection point on the Pareto boundary of the achievable 
rate region with any prescribed "rate profile" vector. This paper shows that such a problem can be transformed into 
a sequence of signal-to-interference-plus-noise ratio (SINR) feasibility problems, which can be solved efficiently by 
existing techniques. Numerical results validate that the proposed algorithms can achieve the global WSR maximum 
for the SISO, SIMO or MISO GIC, which serves as a performance benchmark for existing heuristic algorithms. 

Index Terms 

Beamforming, power control, interference channel, multi-antenna system, non-linear optimization, weighted 
sum-rate maximization. 



I. Introduction 

Gaussian interference channel (GIC) is a basic mathematical model that characterizes many real-life interference- 
limited communication systems. The information-theoretic study on the GIC has a long history, but the capacity 
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region of the GIC still remains unknown in general, even for the two-user case. The best achievable rate region for 
the two-user GIC to date was established by Han and Kobayashi in (J, which utilizes rate splitting at transmitters, 
joint decoding at receivers, and time sharing among codebooks. This achievable rate region was recently proven 
to be within 1-bit of the capacity region of the GIC in |2). However, capacity-approaching techniques in general 
require non-linear multi-user encoding and decoding, which may not be suitable for practical systems. A more 
pragmatic approach that leads to suboptimal achievable rates is to allow only single-user encoding and decoding 
by treating the interference from all other unintended users as additive Gaussian noise. For this approach, the 
key design challenge lies in how to optimally allocate transmit resources such as power, bandwidth, and antenna 
beam among different users to minimize the network performance loss due to their mutual interference. Recently, 
showed that the circularly symmetric complex Gaussian (CSCG) distribution for the transmitted signals is 
in general non-optimal for the rate maximization in GIC with the interference treated as noise. By means of 
symbol extensions over time and/or asymmetric complex signaling, the weighted sum-rate (WSR) of GIC can be 
further improved. However, to our best knowledge, applying such techniques will result in more complicated WSR 
maximization problems, for which how to obtain the globally optimal solutions still remains an open problem, 
even for the case of 2-user GIC. Thus, for simplicity, in this paper we adopt the conventional assumption for the 
GIC that the transmitted signals have an independent CSCG distribution over time. 

The research on the GIC with interference treated as noise has recently drawn significant attention due to the 
advance in cooperative inter-cell interference (ICI) management for cellular networks. Traditionally, most of the 
studies on resource allocation for cellular networks focus on the single-cell setup, while the ICI experienced by 
a receiver in one cell caused by the transmitters in other cells is minimized by means of frequency reuse, which 
avoids the same frequency band to be used by adjacent cells. However, most beyond-3G wireless systems advocate 
to increase the frequency reuse factor and even allow it to be one or so-called "universal frequency reuse", due 
to which the issue of ICI becomes more crucial. Consequently, joint resource allocation across neighboring cells 
becomes a practically appealing approach for managing the ICI. If the mobile stations (MSs) in each cell are 
separated for transmission in frequency via orthogonal frequency-division multiple-access (OFDMA) or in time 
via time-division multiple-access (TDMA), then the active links in different cells transmitting at the same frequency 
tone or in the same time slot will interfere with each other, which can be modeled by a GIC. More specifically, 
if the base stations (BS) and MSs are each equipped with one single antenna, the system can be modeled as the 
single-input single-output (SISO) GIC, termed as SISO-IC. If the BSs are each equipped with multiple antennas 
while MSs are each equipped with one single antenna, then in the uplink the system can be modeled as the single- 
input multiple-output (SIMO) GIC, termed as SIMO-IC, and in the downlink as the multiple-input single-input 
(MISO) GIC, termed as MISO-IcQ 

The achievable rate region of SISO-IC, SIMO-IC or MISO-IC, with the single-user detection (SUD) by treating 
the interference as Gaussian noise, is in general a non-convex set due to the coupled interference among users. As 

'it is worth noting that even for the traditional single-cell setup with space-division multiple-access (SDMA), i.e., the multi-antenna BS 
simultaneously communicating with more than one single-antenna MSs, the MISO-IC and SIMO-IC models are also applicable if the linear 
transmit/receive precoding/equalization is implemented at the BS. 



3 



a result, how to efficiently find the optimal power control and/or beamforming solutions to achieve the maximum 
WSR for different types of GICs is a challenging problem. It is worth noting that a great deal of valuable scholarly 
work [4]-[14] has contributed to resolving this problem. For SISO-IC, various efficient power control schemes have 
been studied. The authors in (H showed that in the two-user case the optimal power allocation to the sum-rate 
maximization problem is "binary", i.e., either one user transmits with full power and the other user shuts down, 
or both users transmit with full power. However, this result does not hold in general when the number of users is 
greater than two. Based on game theory, an "asynchronous distributed pricing (ADP)" algorithm was proposed in 
whereby locally optimal solutions can be obtained for WSR maximization. In |6), the WSR maximization problem 
was transformed into a signomial programming (SP) problem, which was efficiently solved by constructing a series 
of geometric programming (GP) problems through the approach of successive convex approximation. Similar to 
ADP, this algorithm only guarantees locally optimal solutions. As for the case of parallel SISO-IC, the authors 
in 0, 1U showed that the duality gap for the WSR maximization problem is zero when the number of parallel 
GICs becomes asymptotically large. As a result, the Lagrange duality method can be applied to decouple the 
problem into parallel sub-problems in the dual domain. However, the power optimization in each sub-problem for 
a given GIC is still non-convex. For an extensive survey of power control algorithms for SISO-IC, please refer to 
0. Furthermore, for MISO-IC, the optimality of transmit beamforming for achieving the maximum WSR with 
SUD has been proven in ifTOl , ifTTTl . In ifTOl , |[T2l . the complete characterization of all Pareto optimal rates for 
MISO-IC was studied. To maximize the WSR, an iterative algorithm was proposed in lfT3l from an egotistic versus 
altruistic viewpoint, and other "price-based" algorithms (see, e.g., [14] and references therein) were also developed. 
However, these algorithms in general cannot achieve the global WSR maximum for MISO-IC. 

Different from the above prior work in which the power and/or beamforming vectors were optimized directly for 
WSR maximization in the GIC, in this paper we propose a new approach that maximizes the WSR in the achievable 
rate region of the GIC directly. This approach is based on the following two key observations: (1) the WSR is a 
strictly increasing function with respect to individual user rates; and (2) the achievable rate region is a "normal" set 
lTT5l . Accordingly, the WSR maximization problem for the GIC belongs to the class of optimization problems so- 
called monotonic optimization over a normal set, for which the global optimality can be achieved by an iterative 
"outer polyblock approximation" algorithm |fT31 . However, one challenging requirement of this algorithm is a 
unique characterization of the Pareto boundary of the achievable rate region since at each iteration of the algorithm 
one particular point on the Pareto boundary that corresponds to the maximum achievable sum-rate in a prescribed 
direction needs to be determined. This problem is efficiently solved in this paper by utilizing a so-called "rate 
profile" approach [16], which transforms the original problem into a sequence of signal-to-interference -plus-noise 
ratio (SINR) feasibility problems. It is also shown in this paper that such feasibility problems can be efficiently 
solved by existing techniques for various types of GICs. 

It is worth noting that rate profile was first proposed in |[T6l as an alternative method to WSR maximization for 
characterizing the Pareto boundary of the capacity region for the multi-antenna Gaussian multiple-access channel 
(MAC), which is a convex set. This method was later applied to characterize the Pareto boundary of non-convex 
rate regions for the MISO-IC in ifTOTl and the two-way multi-antenna relay channel in ifTTl . for which the WSR 
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maximization approach is not directly applicable. A very similar idea to rate profile was also proposed in 11181 . 
where the proportional rate fairness is imposed as a constraint for WSR maximization in multi-user OFDM systems. 
As for the outer polyblock approximation algorithm, it was first proposed in |fT5l. and later applied in [19 ] and EU1 
to solve the WSR maximization problems for the GIC. In |[T9l , this algorithm was applied for SISO-IC together 
with the generalized linear fractional programming, which, however, cannot be extended to SIMO-IC or MISO-IC. 
In EOl . this algorithm was applied to the two-user MISO-IC by exploiting a prior result in lTT2l that the optimal 
transmit beamforming vector to achieve any Pareto boundary rate -pair can be expressed as a linear combination 
of the zero-forcing (ZF) and maximum-ratio transmission (MRT) beamformers. However, this result only holds 
for the two-user MISO-IC and thus how to extend the algorithm in [20] to MISO-IC with more than two users 
remains unknown. In comparison, in this paper we show that by jointly utilizing the outer polyblock approximation 
algorithm and rate profile approach, the global optimality of the WSR maximization problem can be achieved for 
all SISO-IC, SIMO-IC and MISO-IC, with arbitrary number of users. 

It is also worth noting that for the WSR maximization in SISO-IC, besides |[T9l that applies the outer polyblock 
approximation algorithm, there have been other algorithms developed based on the branch and bound method. 
For example, in Ell and [22], branch and bound methods combined with difference of convex functions (DC) 
programming have been proposed. A generalized branch and bound method applicable to problems in which the 
objective function cannot be expressed in the form of DC, has also been proposed in E3l . In this paper, we propose 
an alternative approach to that in the above prior work, whereby the WSR maximization problems for SISO-IC, 
SIMO-IC, and MISO-IC are all solvable. 

The rest of this paper is organized as follows. Section [II] introduces the system models for various GICs including 
SISO-IC, SIMO-IC and MISO-IC, and formulates their WSR maximization problems. Section In] presents a new 
framework to solve the formulated problems based on monotonic optimization and rate profile techniques. Section 
HVl completes the proposed algorithms by addressing the solutions to various SINR feasibility problems. Section 
|Vl provides numerical examples to validate the proposed results. Finally, Section [VI] concludes the paper. 

Notation: Scalars are denoted by lower-case letters, vectors denoted by bold-face lower-case letters, and matrices 
denoted by bold-face upper-case letters. I and denote an identity matrix and an all-zero matrix, respectively, 
with appropriate dimensions. For a square matrix S, S^ 1 denotes its inverse (if S is full-rank). For a matrix M of 
arbitrary size, M H and M T denote the conjugate transpose and transpose of M, respectively. Diag(Xi, • • • , Xjc) 
denotes a block diagonal matrix with the diagonal matrices given by X\, • • • , Xk- The distribution of a CSCG 
random vector with mean vector x and co variance matrix I] is denoted by CJ\f(x, £); and ~ stands for "distributed 
as". C xxy denotes the space of x x y complex matrices. R denotes the real number space, while M. x denotes the 
x x 1 real vector space and R+ denotes its non-negative orthants. ||a;|| denotes the Euclidean norm of a complex 
vector x. denotes a vector with its kth component being 1, and all other components being 0. For two real 
vectors x and y, x > y means that x is greater than or equal to y in a component-wise manner. 

II. System Model and Problem Formulation 

We consider a K-user GIC, in which K mutually interfering wireless links communicate simultaneously over 
a common bandwidth, as shown in Fig. [T] Firstly, consider the case where all transmitters and receivers are each 
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Fig. 1. System model for the K-usei MISO-IC (SISO-IC if each transmitter has one single antenna, or SIMO-IC in the reverse link 
transmission). 

equipped with one single antenna. The system is thus modeled as SISO-IC, for which the discrete-time baseband 
signal received at the Mi receiver is given by 

Vk = h k ,kVPk~ x k + hk dVPJ x J + z ki k = K, (1) 

where hkj denotes the complex channel gain from the jth transmitter to the kth receiver, pk denotes the transmit 
power of the Mi transmitter, Xk denotes the transmitted signal from the kth transmitter, and Zk denotes the 
background noise at the kth receiver. It is assumed that Zk ~ CA/"(0, cr^), V/c, and all z^'s are independent. 

We assume independent encoding across different transmitters and thus x^'s are independent over k. It is also 
assumed that the Gaussian codebook is used and thus Xk ~ CAf(0, 1). Accordingly, the SINR of the Mi receiver 
is expressed as 

SISO-IC _ \\hk,k\\ 2 Pk n , 

lk " £IIMI 2 ^ + ^' U 

Remark 2.1: It is worth noting that in the above signal model, we have made the following two assumptions: 
Al. The interference is treated as additive Gaussian noise. 

A2. The Gaussian input xu for user k is assumed to be CSCG distributed and independent over time, i.e., 

asymmetric Gaussian signalling with time-domain symbol expansion in is not used. 
Note that for the subsequent studies on SIMO-IC and MISO-IC, the above two assumptions are similarly made. 

Secondly, consider the case where all transmitters are each equipped with one single antenna but each receiver is 
equipped with multiple antennas, i.e., SIMO-IC. Assuming that the Mi receiver is equipped with Mk > 1 antennas, 
its discrete-time baseband received signal is given by 

yk = w^(h kj kVPkXk + ^2hkj^/p]xj + Zk), k = l,---,K, (3) 
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where G C lxMfc is the receive beamforming vector for the fcth receiver, h k j G C MfeXl is the channel vector 
from the jth transmitter to the fcth receiver, and z k G <c MfcXl is the noise vector at the kth receiver. It is assumed 
that z k ~ CM(0, ^1). Thus, the SINR of the kth receiver can be expressed as 

1 1 H L. 1 1 2 

simo-ic Vk\\w k n k:k \\ 

Thirdly, consider the MISO-IC case in which all transmitters are each equipped with multiple antennas while 
each receiver is equipped with one single antenna. Assume that the kth transmitter is equipped with N k > 1 
antennas. The discrete-time baseband signal at the kth receiver is then given by 

Vk = h% ik v k x k + ^2 h k,j v J x j + z k, k = 1, ■ ■ ■ , K, (5) 

where v k G <C NkXl is the transmit beamforming vector at the kth transmitter, and h k j G C lxNj denotes the 
channel vector from the jth transmitter to the kth receiver. Accordingly, the SINR of the kth receiver can be 
expressed as 

\\uH [ 1 2 
MISO-IC ll^fc.fc^fcll 

jj^k 

Remark 2.2: It is worth noting that for the above signal model for MISO-IC, we assume that all transmitters 
employ rank-one beamforming. This is because it has been shown in ifTOTl and iTTTTl that under Assumptions Al 
and A2, beamforming achieves all the points on the Pareto boundary of the achievable rate region for MISO-IC, 
i.e., beamforming is optimal for WSR maximization in MISO-IC. 

With 7fc defined in ©, (0]) or ((6]), the achievable rate of the kth receiver can be formulated as 

flfc(7fc)=log 2 (l + 7fc)i k = !,■■■ ,K. (7) 
Next, we define the achievable rate region for each type of GIC, which constitutes all the rate-tuples simultane- 
ously achievable by all the users under a given set of transmit-power constraints denoted by p™^, . . . , P™ ax : 

^ SIS °- IC 4 (J [(n, . . . ,r K ) : < r k < R k (^°~ lC ),k = 1,...,k\, (8) 



^ SIM °- IC ^ U { in> ■ ■ ■ > *k) ■ < n < i? fc (7f IM °- IC ), k = 1, . . - , A , (9) 

{p k },{W k }:0<p k <P™^, Vfc ^ > 

^ MIS °- IC 4 J . . . ,r K ) : < r k < R k ( 7 ^%k = I,..., A. (10) 

{V h }:0<\\V k \\*<PF**, Vfe ^ ' 

The upper-right boundary of each defined rate region is called the Pareto boundary, constituted by rate-tuples for 
each of which it is impossible to improve one particular user's rate without decreasing the rate of at least one of 
the other users. 

The WSR maximization problems for SISO-IC, SIMO-IC and MISO-IC are then formulated as (P1.1)-(P1.3) as 
follows. 

K 

(Pl.l): Maximize U{p) := X^fcflfcfr? S °~ IC ) 
P k=i 

Subject to < p k < P fc max , Vk, 
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K 

(P1.2) : Maximize U(W,p) := V ^Rkilf MO ~ lc ) 
w P iTi 

Subject to < p k < P fc max , V/c, 

K 

(P1.3) : Maximize U{V) := J] » k Rk{lf l *°~ lQ ) 
V k=i 

Subject to \\v k \\ 2 < P fc max ,VA;, 

where p = (p%, ■ ■ ■ ,Pk) denotes the transmit power vector, W = (wi, • • • , wk) and V = (v\, • • • , vk) constitute 
the receive and transmit beamforming vectors, respectively, and p,^ is the non-negative rate weight for user k. Since 
the objective functions are all non-concave with respect to the power values or beamforming vectors due to the 
coupled interference, all the WSR maximization problems in (P1.1)-(P1.3) are non-convex and thus cannot be 
solved globally optimally by conventional convex optimization techniques. 

III. Proposed Solution Based on Outer Polyblock Approximation and Rate Profile 

In this section, we solve the formulated WSR maximization problems in (P1.1)-(P1.3) globally optimally by a 
new approach based on the outer polyblock approximation and rate profile techniques. 

A. A New Look at the Problem: Optimizing WSR Directly in Rate Region 

Traditionally, Problems (P1.1)-(P1.3) are solved in the power allocation and/or beamforming domain, which 
results in non-convex optimization problems. In this subsection, we study the WSR maximization problem utilizing 
a new formulation, which maximizes the WSR directly in the achievable rate region. 

If the achievable rate vector r = (R\, • • • , Rk) is treated as the design variable, where is the achievable 
rate of user k defined in (O, the WSR maximization problems (P1.1)-(P1.3) can be unified in the following form. 

K 

(P2) : Maximize U(r) := ^//fc-Rfc 

fc=i 

Subject to r 6 1Z, 

where the rate region TZ is defined in ®, © or flO]) for SISO-IC, SIMO-IC or MISO-IC. 

Next, we will show that Problem (P2) belongs to one special class of optimization problems: monotonic 
optimization over a "normal" set. Two useful definitions are given first as follows. 

Definition 3.1: A function / : M. n — > R is said to be strictly increasing on WL if for any x',x G R", x' > x 
and x' ytz x imply that f(x') > f(x). 

Definition 3.2: A set V G W]_ is called normal if given any point x G V, all the points x' with < x' < x 
satisfy that x' G V. 

Based on the above definitions, we declare the following two facts regarding Problem (P2), which can be easily 
verified to be true. 

Fact 1: The objective function of Problem (P2) is a strictly increasing function with respect to r. 
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Fact 2: The achievable rate region defined in ([U), © or (fTOl) is a normal set. 

Facts [T] and |2] imply that Problem (P2) maximizes a strictly increasing function over a normal set. In 031. the 
"outer polyblock approximation" algorithm was proposed to achieve the global optimality for this type of problems. 
In the following, we will apply this algorithm to solve Problem (P2). 

B. Outer Polyblock Approximation Algorithm 

In this subsection, we introduce the outer polyblock approximation algorithm to solve Problem (P2). First, two 
definitions are given as follows. 

Definition 3.3: Given any vector v £ R™, the hyper rectangle [0,v] = {x\0 < x < v} is referred to as a box 
with vertex v. 

Definition 3.4: A set is called a polyblock if it is the union of a finite number of boxes. 
Next, we show one important property of the polyblock in the following proposition. 

Proposition 3.1: The maximum of a strictly increasing function f(x) over a polyblock is achieved at one of 
the vertices of the polyblock. 

Proof: Suppose that x* is the globally optimal solution over the polyblock, and it is not a vertex. Then, 
there exists at least one vertex x' satisfying x' > x* but x' ^ x*. Since f(x) is a strictly increasing function, 
f(x*) < f(x') must hold, which contradicts to the presumption. The proof is thus completed. ■ 

According to Proposition [3J] the maximum of an increasing function over a polyblock can be obtained efficiently 
by enumeration of the vertices of that polyblock. Consequently, we can construct a sequence of polyblocks to 
approximate the rate region 1Z with the increasing accuracy for Problem (P2). In other words, we need to find an 
iterative method to generate a sequence of polyblocks of shrinking sizes such that 

pW D p(2) 3 ... D n, (ii) 

lim [ max U(r)] = maxf/Yr), (12) 
n-s.oo L reP<™> veil 

where p( n ) denotes the polyblock generated at the nth iteration. 

Next, we present one method to generate the polyblocks satisfying (fTTb and (TT2l) . Let denote the set 

containing all the vertices of the polyblock P^ n \ The vertex that achieves the maximum WSR in polyblock P^ 

can be formulated as 

~ z (n) = arg max TJ( Z \ ( 13 ) 

Define 5z^ as the line that connects the two points and zS n \ and as the intersection point on the Pareto 
boundary with the line dz^. The following method can be used to generate K new vertices adjacent to z^: 

z W.* = 5(»)_( 5 W_ r W) Cf) i = l,...,K, (14) 

where z( n )' J denotes the ith new vertex generated at the nth iteration; and r„- denote the ith element of 
vectors zS n ^ and r^ n \ respectively. Then, the new vertex set can be expressed as 

Z (n+1) = Z {n)y z {n) \J {z (n),^ . . . ; z (n),Ky (15) 
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The polyblock P <n) which R, 
approaches the rate region 



with increasing n 



Fig. 2. Illustration of the procedure for constructing new polyblocks. 



Each vertex in the set z( n+1 ) defines a box, and thus the new polyblock is the union of all these boxes. 

An illustration about the above procedure to generate polyblocks for the case of two-user rate region is given in 
Fig. [2] In the following proposition, we show the feasibility of the above polyblock generation method. 

Proposition 3.2: If the rate region 1Z is a normal set (as we have already shown), the polyblocks generated by 
(fT3T ) satisfy (fTTT ). 

Proof: Please refer to ITT51 . ■ 
Proposition 13.21 ensures that the above polyblock generation method can be used to approximate the rate region 
from the outside with increasing accuracy. Let r* = (R*, • • • , R* K ) denote the optimal solution to Problem (P2). 
Based on the above method, in the following we present an algorithm to find r* in the rate region 1Z. It is worth 
noting that r* must be on the Pareto boundary of the rate region; thus, we only need to search over the Pareto 
boundary to find r*. 

The outer polyblock approximation algorithm works iteratively as follows. In the nth iteration, the optimal vertex 
zS n ^ is first obtained by ( fl"3l) . According to Proposition 13. II in the polyblock p( n ) the maximum WSR is U(z^). 
Since Proposition 13.21 implies that always contains the rate region is an upper bound of U{r*). 

Then, the intersection point r^ n ' on the Pareto boundary with the line dz^ is obtained. Define the best intersection 
point up to the nth iteration as 

f (n) = argmax{[/(7»), U{r {n ~ l) )}. (16) 
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Consequently, U(f^) is the tightest lower bound of U(r*) by the nth iteration. Next, we can compute the value 
of U(z^) — U(r^), which is the difference between the upper and lower bounds of the optimal value of Problem 
(P2) at the nth iteration. If this difference is less than r\ (a small positive number), the algorithm can terminate 
and f ( n ' is at least an //-optimal solution to Problem (P2) because 

U(r*) - U(r {n) ) < U(z^) - U(r {n) ) < rj. (17) 

Otherwise, we construct a new polyblock p( n+1 ) by the above polyblock generation method. We repeat the above 
procedure until an ^-optimal solution is found. 

TABLE I 

Algorithm!) Outer Polyblock Approximation Algorithm for Solving Problem (P2) 

a) Initialize: Set n = 1, Z (1) = {z (1) }; 

b) While (e, 77) -accuracy is not reached, do 

1) Find the optimal vertex zS n ^ that maximizes the WSR in the set Z^ based on 

zr n ' — arg max U(z), (18) 

where e is a small positive number and Z^ = {z £ Z^\z^ > e, VA;}; 

2) Compute the intersection point r'"' on the Pareto boundary of the rate region 1Z with the line Sz^ ; 

3) Update the best intersection point until the nth iteration f'™' according to 0161 ); 

4) If U(z (n) ) - U(r (n) ) < r), then 

Stop and f'"' is an (e, r7)-optimal solution to Problem (P2); 

5) else 

Compute K new vertices that are adjacent to zS n ' by d 14b and update the vertex set Z^ n+1 "> by j 1 5b : 

6) end 

7) n = n + 1; 

c) end 



The above algorithm, denoted as Algorithm H is summarized in Table U It is worth noting that in Algorithm 
H is obtained by enumeration in the set 2^ n ' rather than Z^ n \ This is because in |fT31 it was shown that if 
the optimal solution lies in a strip defined by {r*|0 < R* k < e} with arbitrary k and a small value e > 0, then as 
zS n ^ approaches this strip, the algorithm converges very slowly. Consequently, e is chosen to balance the tradeoff 
between the accuracy and complexity of Algorithm U With e, Algorithm U solves the following problem 

K 

(P2 - A) : Maximize U(r) := ^/j, k R k 

k=l 

Subject to r£ TZ e , 

where Vf is defined as 

K e = Kn{{r x ,--- ,r K ) :r k >e,Vk}. (19) 
Thus, the corresponding solution is called an (e, ?])-optimal solution to Problem P2. 
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Next, we address the convergence issue of Algorithm J] According to Proposition 13.21 P^ D p( n + l ) always 
holds. Moreover, the optimal vertex z^> is removed from after each iteration. Thus, U(z^ n+1 ^) < U(z^) 

also holds. Furthermore, the lower bound U(r^) is non-decreasing. Consequently, the value of U(z^) — U(r^) 
will decrease after each iteration. It was shown in lff5ll that as n increases, the difference between the upper and 
lower bounds can be reduced to an arbitrary small positive number in a finite number of iterations. Thus, Algorithm 
U converges given small positive values e and n. More details about the selection of the values of e and n will be 
given later in Section |V-Aj 

Last, we explain how to obtain an initial vertex z^ = (z^,--- ,z^)) for the first iteration of Algorithm J] 
Since the box [0, «W] needs to contain the rate region 1Z, for any user k, zip can be obtained when all other users 
switch off their transmission (thus no interference exists for user k), and user k transmits its maximum power 
pmax More specifically, for SISO-IC, 

pmax II/, 1 1 2 

4 1} = log 2 (l + ^ '7'*" ), V*. (20) 



Since for MISO-IC, 



MISO-IC _ \\hk,k w k\\ 2 ll^fcjg^fcll 2 

< ii fc n ii Ml! < _fe — ii *.*» V fe, (21) 



where (a) is due to the Cauchy-Schwarz inequality, zP can thus be set as 



pmax \\u 1 1 2 

g> = log 2 (l + fc "**»*" ), Vfc. (22) 



The initial vertex for SIMO-IC can be obtained similarly to (1221 . and is thus omitted for brevity. 
To summarize, the only challenge that remains unaddressed in Algorithm U is on how to compute the intersection 
point r( n ) on the Pareto rate boundary with the line Sz^ at the nth iteration, which will be addressed next. 

C. Finding Intersection Points by the "Rate Profile " Approach 

In this subsection, we show how to obtain the intersection point on the Pareto boundary of the rate region with 

K K 

the fine 5z^ n \ to complete Algorithm HI Let i? sum = ^2 Rk denote the sum-rate of all the users, a = z {n) / ]T z^ n) 

k=l fc=1 

denote the slope of the line 5z^ n \ Consequently, the intersection point at the nth iteration can be expressed as 
r ( n ) = R* um a, where R* nm is the optimal value of the following problem: 

Maximize i? sum 
Subject to -Rsum 

at £ 1Z. (23) 

The above approach to find the intersection point on the Pareto boundary of the rate region is known as rate 
profile ifTOl , lfT6l , ifT/Tl . In the following, we solve Problem (|23l to obtain the intersection point on the Pareto 
boundary with a given 5z^ n \ 
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Problem (1231 ) is solvable via solving a sequence of feasibility problems shown as follows. Given a target sum-rate 
-Rsum, the feasibility problems for SISO-IC, SIMO-IC and MISO-IC can be expressed in the following problems 
(P3.1)-(P3.3), respectively. 

(P3.1) : Find {p k } 

Subject to log 2 (l +7fc IS °~ IC ) > a k R sum , V/c 
Pk<Pr x , Vfc. 

(P3.2) : Find {w k }, { Pk } 

Subject to log 2 (l + 7fc IM °" IC ) > a k R sum , V/c 

(P3.3) : Find {v k } 

Subject to log 2 (l + 7 f ISO ~ IC ) > a k R sum , V/c 

||^||2 < pn« j 

If any of Problems (P3.1), (P3.2) and (P3.3) is feasible, it follows that i?g Um > -Rsum; otherwise, i?s Um < ^sum- 
Hence, i?sum can be obtained for Problem (1231 by applying a simple bisection method l24l . for which the detail 
is omitted for brevity. 

The remaining challenge is on solving the feasibility problems (P3.1)-(P3.3), which is addressed next. Let 
% = 2 akRsura — 1, V/c. Then, the first constraint of each feasibility problem can be re-expressed as 

7fc>7fe, V*. (24) 

Therefore, given any sum-rate target R sum , the feasibility problems (P3.1)-(P3.3) are equivalent to finding whether 
a corresponding SINR target vector 7 = (71, • • • ,jk) is achievable. In the next section, we will propose efficient 
algorithms to solve these SINR feasibility problems. 

Remark 3.1: In the case where a set of minimum rate constraints R k > i?™ in , V/c, are added to the WSR 
maximization problem (P2), where i?™ in is the minimum rate required for user k, we can solve this new problem 
by modifying Algorithm U as follows. Since the new rate region 1Z' is the intersection of the original rate region 
with the set {(r%, • • • ,tjc) '■ r k > ^/J 1111 ; V/c}, we should change the initial point from to r min in Algorithm 
H where r mm = (i? mm , • • • , i? mm ) is the rate constraint vector. Thus, at each iteration we need to compute the 
intersection point on the Pareto boundary with the line passing through the optimal vertex and the point r min 
(instead of in Algorithm |I]). In addition, any point r on this line with R sum = Y2k=l ^-k can be rewritten as 

K 

r = r min + a(i? sum - £ Rf»), (25) 

k=l 

where the rate profile a is obtained by a = ^'J — — . 
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IV. Solutions to SINR Feasibility Problems 

In this section, we solve Problems (P3.1)-(P3.3) subject to the equivalent SINR constraints given in (l24b for 
SISO-IC, SIMO-IC and MISO-IC, respectively. 

A. The SISO-IC Case 

We first study the feasibility problem (P3.1) for SISO-IC. Given a SINR target vector 7 = (71, • • • ,jk) with 
7^ = 2 akRsuni — 1, we will check whether it is achievable under users' individual power constraints. 
Let G denote the K x K normalized channel gain matrix given by 

G J IMP ' ^ (26) 
[ 0, fc = j, 

and 77 denote the K x 1 normalized noise vector given by 

To achieve the SINR target, the transmit power vector for users is given by 

p=(I-G) 1 rj. (28) 

Let p(B) denote the spectral radius (defined as the maximum eigenvalue in absolute value) of the non-negative 
matrix B. The following propositions were shown in |[25l . which play important roles in solving Problem (P3.1). 

Proposition 4.1: The power allocation p given by (l28l) satisfies p > if and only if p(G) < 1. 

Proposition 4.2: If p(G) < 1, the power allocation p given by (l28l) is component-wise minimum in the sense 
that any other power allocation p' that satisfies (1241 ) needs to satisfy p' > p. 

Propositions 14. 11 and 14.21 imply that a SINR target vector 7 is feasible if and only if: (1) p(G) < 1, and (2) the 
power solution obtained by (l28l) satisfies p^ < P™ ax , V/c. Consequently, we propose Algorithm [TT] in Table [TT] to 
solve Problem (P3.1). 

TABLE II 

AlgorithmJII Algorithm for Solving Problem (P3.1) 

a) Given any SINR target vector 7 = (71, • • • , 7jf), compute the spectrum radius of matrix G. If it is larger than 1, conclude that there 
is no feasible power allocation to meet the SINR target and exit the algorithm; otherwise, go to step b); 

b) Compute the power allocation p by 128t . and check for any user k, whether the power constraint pt < P™ ax is satisfied. If so, 
conclude that the SINR target is feasible; otherwise, the SINR target is not feasible. 



Remark 4.1: It is worth comparing Algorithm U with the MAPEL algorithm proposed in [19]. MAPEL solves 
Problem (Pl.l) for SISO-IC in the SINR region (as opposed to the rate region in our approach) due to the fact 
that the problem to characterize the Pareto boundary of the SINR region for SISO-IC can be transformed into a 
generalized linear fractional programming problem and thus efficiently solved by Dinkelbach-type algorithm j26l . 
However, this transformation does not work for SIMO-IC or MISO-IC if the beamforming vectors are involved. 
Consequently, MAPEL cannot be extended to the GIC with multiple antennas. As comparison, in this paper we 
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solve the WSR maximization problem in the rate region directly because the Pareto boundary can be characterized 
completely by the rate profile approach, as along as the associated SINR feasibility problem can be solved. Thus, 
our proposed algorithm is more applicable than MAPEL in solving the WSR maximization problems for SIMO-IC 
and MISO-IC, as shown next. 

B. The SIMO-IC Case 

The feasibility of Problem (P3.2) can be checked by using the optimal value of the following SINR balancing 
problem: 

Ik 

Maximize min — 

l<k<K 7fc 

Subject to Pk < P fc max , VJfe. (29) 

If the optimal value of Problem (1291 is no smaller than 1, then the SINR target vector 7 = (71,- •• ,7if) is 
achievable; otherwise, this SINR target cannot be achieved. 

In (271, an efficient algorithm was proposed to solve a SINR balancing problem similar to Problem d29l ), where 
only one sum-power constraint is imposed. However, the algorithm in 1271 cannot be directly applied to solve 
Problem d29l ) due to multiple users' individual power constraints. To utilize the algorithm proposed in ||27l , we 
decouple Problem d29l into K sub-problems, with the ith sub-problem formulated as: 

Ik 

Maximize min — 

l<k<K 7^ 

Subject to K < i^ max . (30) 

Therefore, for the ith sub-problem only the ith user's power constraint is considered. Next, we show how to solve 
Problem 001 by extending the algorithm in 12711, and then reveal an important relationship between Problems d29l 
and 001 . based upon which we further propose an efficient algorithm to solve Problem d29l . 
I) Solution to Problem A301) : 

In this part, we extend the algorithm proposed in |[27l to solve Problem 001 for a given i. 
One important property of the SINR balancing problem in 001 is that given any receive beamforming vectors 
W = (wi, ■ ■ ■ ,wk), the corresponding optimal power allocation p must satisfy the following two conditions: 

lk(W, P ) =g(T y ) vfc (31) 

Pi = P t max , (32) 

where C(W) is the maximum SINR balancing value for all users given W. 

We justify the above conditions as follows. (OTI ) can be shown by contradiction. Supposing that the SINR 
balancing values are not the same for all the users, then we select the user with the highest SINR balancing value 
and decrease its transmit power by a small amount such that its new SINR balancing value is still above min 

k ^ k 

Since the other users' SINR balancing values will increase, the minimum SINR balancing value among all the 
users will increase accordingly. Thus, whenever the SINR balancing values are not the same for all users, we can 
proceed as above to improve the optimal value. Hence, 0TT ) must hold. Similarly to show 02b by contradiction, 
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pmax 

suppose pi < Pj max . With a = > 1, we can multiply the transmit power values of each user by a, and the 



SINRs of all users will be increased accordingly. Hence, (1321 ) must hold. 
We can express (|3TT) for all fc's in the following matrix form: 

p— L-=D^(W)p + Da, (33) 

where D = Diag{ — — -, • • • , H T — r^}, <x = [aj ||"»I>i|| 2 , • • • , a 2 K \\w K \\ 2 ] T , and the K X K non-negative 

II W 1 »tl,l|| \\W k IIk,K\\ 

matrix *&(W) is a function of W defined as 



Mw)] kJ 



\\ w k h k,j\\ 2 , k^j 
0, k = j. 

By multiplying both sides of (l33l by ef, we obtain 



(34) 



j pmax 

^ciWi^c^)^- *^^'^' <35) 



Therefore, by combining (1331) and (1351) . it follows that 

1 



C(W) 



Pext = Ai(W)p cxt , (36) 



where p ext = | P | and 




A(^) = 1 . (37) 



Next, we show one important property for (1361 ) in the following lemma. 

Lemma 4.1: Given any fixed W, there exists a unique solution (p, C(W)) to the equation in 

Proof: Please refer to Appendix El ■ 

According to Perron-Frobenius theory [29], for any nonnegative matrix, there is at least one positive eigenvalue 
and the spectral radius of the matrix is equal to the largest positive eigenvalue. Furthermore, according to Lemma 
14.11 there is only one strictly positive eigenvalue to matrix Aj(W'). Accordingly, it follows from (l36l ) that given 
W, the inverse of the optimal SINR balancing value 1/C(W) is the spectral radius of Ai(W). Consequently, the 
maximum SINR balancing solution to Problem (l30l ) is obtained as 



Next, by defining a cost function as 



° min p(Ai(W)Y (38) 
W 



T(W, Pcxt ) = max^^. (39) 
x>o x T p ext 

then the min-max characterization of the spectral radius of Aj(W) can be expressed as E71 . ll30l 

p(MW)) = mmT(W, Pext ). (40) 

■I cxt 



Taking (gOj) into (|38ll . it follows that 

1 



minminT(W,p ext ). (41) 
C* W P oxt 
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Similar to E71 . we can solve Problem (|4TT ) via the alternating optimization shown as follows. First, given W, 
we find the optimal power allocation for p ext . Let p ext denote the dominant eigenvector corresponding to the 
spectral radius of Ai(W). It then follows that 



XTA ^T PeXt = P(MW)) = mmT(W,p cxt ). (42) 

x Pext f'ext 

Thus, p ext is the optimal power allocation given W. 

Furthermore, we know that given any power allocation p ext , the optimal receive beamformer in W to maximize 
the SINR is minimum-mean-squared-error (MMSE) based for each of the users. Therefore, we propose an iterative 
algorithm in Table [TTll denoted as Algorithm JIIJ to solve Problem (|30l ). 

TABLE III 

Algorithm IIIII Algorithm for Solving Problem ( T30T > 

a) Initialize: n = 0, p (0) = [0, • • • , 0] T and p (0) = oo; 

b) repeat 

1) n = n + 1; 

2) Update W {n) by w[ n) = (£ pj n_1) h fejj h^ + ogl)- 1 !*,*, Vfc; 

3) Update p£™' as the dominant eigenvector of the matrix Ai(W < - n ^); 

4) = P {A r {W (n) )) and C (n) = ^iy; 

c) until p'"- 1 ' - p (n) < e. 



The convergence of Algprithm HIT] can be shown in the following way. Since given any power allocation p^ 
for the nth iteration, W^( n+1 ) minimizes T(W ,p^ t ), i.e., 

T(W^ n+1 \p { :i) < T(WW, P W) = pH. (43) 

Moreover, given V^ n+1 \ p^ 1 minimizes T(VF^ n+1 \p cxt ) as 

p {n+l) = T(W (n+l) )p («+l) ) < T(W(" +1 ),pS). (44) 

Hence, we can guarantee p( n+1 ) < p( n ) after each iteration. Since p is lower-bounded by 0, Algorithm HIT] thus 
converges. 

Finally, the convergence of Algorithm JIlJ to the global optimality of Problem ( f3Qb can be proven similarly 
as Section IV. A in ll27l . and the proof is thus omitted for brevity. After convergence, *y is the maximum 
achievable SINR vector and p( n \ are the optimal power and receive beamforming vectors to achieve this 

SINR vector, respectively. 

2) Solution to Problem (I29D : 

Next, we show that Problem (|29l ) can be efficiently solved via solving Problem (l30l ) for all z's. Let W* and 
p* denote the optimal beamforming vectors and power allocation for Problem d29l ), respectively. Let W* and p* 
denote the optimal beamforming vectors and power allocation for the ith sub-problem in (l30b . respectively. Next, 
we provide a theorem to reveal the relationship between the optimal solutions to Problems (1291 ) and (|30l ). 
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Theorem 4.1: For all sub-problems in (l30l with i = 1, • • • ,K, there exists one and only one sub-problem for 
which the optimal power solution satisfies all users' individual power constraints of Problem (|29l ). Furthermore, 
let i* denote the index of the corresponding sub-problem in 001 ), then it holds that W* = W*,, and p* = p*„. 
Proof: Please refer to Section IV.B in IlTTI . ■ 

Theorem 14. 1 1 reveals that Problem d29l ) can be solved as follows. First, we apply Algorithm UlTI to solve Problem 
(l30l in the order of i = 1, • • • , K. If the optimal power solution to any of these problems satisfies all users' 
individual power constraints, the algorithm terminates, and the obtained optimal power and beamforming solutions 
to Problem d30l are also those to Problem ( f29l ). The above algorithm, denoted by Algorithm |IVJ is summarized in 
Table [IVJ 

TABLE IV 

AlgorithmJvJ Algorithm for Solving Problem J29T ) 



a) Initialize: i = 0; 

b) repeat 

1) i = i + l; 

2) Solve the ith sub-problem in i30\ by Algorithm [HI] and find the optimal beamforming solution W* and power solution p*\ 

3) Check whether p* satisfies all power constraints of Problem l |29t . If so, exit the algorithm and set W* and p* as the optimal 
solution to Problem j29\ ; otherwise, continue the algorithm; 

c) until i — K. 



C. The MISO-IC Case 



In this subsection, we show how to solve the feasibility problem (P3.3) for MISO-IC under the equivalent SINR 
constraints given by (T24)) . It was shown in [10] that this problem can be transformed into a second-order cone 
programming (SOCP) problem, which is briefly described as follows for the sake of completeness. The SINR 
constraints in Problem (P3.3) can be rewritten as 

1 K 

(l + ^)\\h% k v k \\ 2 >^\\h% jV:i \\ 2 + al Vfc. (45) 



Ik 



Without loss of generality, we can assume that h kk Wk is a positive real number, \fk. Thus we can reformulate 
the above SINR constraints as 



1 



1 + ^h% k v k > 

Ik 



K 



(46) 



Denote x = [v 



1 ! 



vj(, 0] of dimension {K 2 + 1) x 1, = [0, ■ ■ ■ ,0, cr fc ] r of dimension (K + 1) x 1, and 



El 



Diag(h kl , ■ ■ ■ , h kK , 0) of dimension (K + 1) x (K 2 + 1), VA;. We further define as 



\KxK 



n KxK jKxK r\KxK 



,0 KxK ,0 Kxl 



(47) 
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where o KxK and Kxl denote the K x K all-zero matrix and K x 1 all-zero vector, respectively, and J KxK 
denotes the K x K identity matrix. Thus, d46l ) can be reformulated as 



\E k x + n k \\ < Jl + —hU k L k x, Vfc. (48) 
V Ik ' 



Moreover, we can reformulate the power constraints as 



\\L k x\\ < VPT^i Vfc - ( 49 ) 

Using (|48T ) and (|49l . Problem (P3.3) can be transformed into a SOCP feasibility problem over x and efficiently 
solvable by existing software IT321 . 

Remark 4.2: It is worth comparing our proposed algorithm with that in EDI for solving the WSR maximization 
problem (PI. 3) for MISO-IC. The algorithm in [20] is based on a prior result in lfT2l that for the special case of 
two-user MISO-IC, any point on the Pareto boundary of the rate region can be achieved by transmit beamforming 
vectors that are obtained by linearly combining the ZF and MRT beamformers. In EDI , the outer polyblock 
approximation algorithm was applied to find the optimal beamformer combining coefficients. However, since this 
result does not hold for MISO-IC with more than two users, the algorithm in E01 cannot be extended to the 
general K-user MISO-IC with K > 2. In contrast, our proposed algorithm can be applied to MISO-IC with 
arbitrary number of users. 

V. Numerical Results 

In this section, we provide numerical results to validate the proposed algorithms in this paper. We assume that 
ii k = 1, V/c, i.e., the sum-rate maximization problem is considered. We also assume that p™ ax = 3, For 
SIMO-IC and MISO-IC, we further assume that M k = 2 and N k = 2, respectively, \/k. The numerical results with 
related discussions are presented in the following subsections. 

A. Convergence Performance 

Firstly, we study the convergence performance of Algorithm U for SISO-IC. We assume that there are 4 users, 
i.e., K = 4, and there is a minimum rate constraint for each user with i?™ m = 0.5, \/k. We set the parameters to 
control the accuracy of Algorithm U as e = 0.01 and 77 = 0.5. We consider the following matrix: 

0.4310 0.0022 0.0105 0.0042 

0.0200 0.4102 0.0180 0.0035 

0.0210 0.0200 0.5162 0.0112 

0.0210 0.0021 0.0063 0.3634 

with each element denoting the power of the corresponding channel gain, i.e., H k j = \\h k j\\ 2 . 

Fig. [3] shows the convergence of Algorithm U under the above channel setup. It is observed that this algorithm 
takes about 300 iterations to converge. The converged sum-rate is 11.4605 with users' individual rates given by 
[3.1982, 2.6297, 2.8441, 2.7884]. To verify that the global sum-rate maximum is achieved, we compare the obtained 
maximum sum-rate with that by an exhaustive search, which is equal to 11.5349. Thus, Algorithm U does achieve 



H 



(50) 
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Fig. 3. Convergence performance of Algorithm Q] for SISO-IC with weak interference channel gains. 




Fig. 4. Convergence performance of Algorithm Q] for SISO-IC with strong interference channel gains. 
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the global optimality of sum-rate maximization within a guaranteed error 11.5349 — 11.4605 = 0.0744, which is 
smaller than the set threshold r\ = 0.5. 

Next, we consider a SISO-IC with stronger cross-user interference channel gains than those given in d50l ) by 
keeping all diagonal elements of H unchanged, but scaling all off-diagonal elements by 10 times. As shown in 
Fig- SI for this new channel setup, Algorithm U takes about 2900 iterations to converge. The converged sum-rate in 
this case is 5.1184 with users' individual rates given by [0.5408, 1.9119, 0.5060, 2.1597], while that obtained by the 
exhaustive search is 5.1392. Thus, as compared to the previous case with weaker interference channel gains, the 
global optimality for sum-rate maximization is achieved in this case with a much slower convergence. The reason 
is as follows. With stronger interference channel gains, the optimal power allocation for sum-rate maximization is 
more likely to render some users transmit at their minimum required rates (e.g., user 1 and user 3 in this example). 
Hence, the corresponding optimal rate values will lie in the strip defined by {r*\R^ m < < R™ m + e} for 
some fc's. Since in Algorithm U each new polyblock is generated from the previous one by cutting off some 
unfit portions, the cuts become shallower and shallower as zS n ^ approaches the above strip. This can be observed 
from Fig. [4] that after the 1300th iteration, the best intersection point has never changed. However, to make 
U(z^) — U{r( n ') < rj hold, another 1700 iterations are taken just to reduce the value of U{z^ n '). Since this 
reduction becomes very inefficient near the strip, the algorithm converges much more slowly to the desired accuracy 
with the increasing of interference channel gains. From this observation, we infer that the values of e and rj need 
to be properly set to balance between the accuracy and convergence speed of our proposed algorithm. 

Next, we give another example to illustrate the important role of parameter e in balancing between the accuracy 
and convergence speed of our proposed algorithm. We assume that K = 3, and there are no minimum rate 
requirements for the users. We consider the following channel matrix: 

0.4310 0.0187 0.0893 " 
H 0.1700 0.4102 0.1530 , (51) 
0.1785 0.1700 0.5162 

with Hfcj = \\hkj\\ 2 . By an exhaustive search, the optimal sum-rate is obtained as 4.8079 with users' individual 
rates given by [3.2146, 1.5933,0]. 

Table [V] shows the convergence speed and the converged sum-rate of our proposed algorithm for different values 
of e with r\ = 0.2. We observe that as e increases, the algorithm convergence speed improves rapidly, but the 
converged sum-rate decreases. When e = 0.45, the difference between the optimal sum-rate and converged sum- 
rate is 4.8079 — 4.5880 = 0.2199, which is even larger than r\ = 0.2. This is because that as we show in Section 
IIII-B I with non-zero e, we are in fact solving Problem P2-A instead of the original problem P2. Consequently, the 
proposed algorithm can only guarantee that the difference between the maximum sum-rate of Problem P2-A and 
the converged sum-rate is less than rj, but not necessarily for Problem P2. Thus, if the value of e is selected to be 
too large such that all the ^-optimal solutions lie in the excluded strips, the difference of the converged sum-rate 
and the maximum sum-rate of Problem (P2) will be larger than rj. Therefore, the value of e should be carefully 
selected based on the value of rj. In this numerical example, we can select e = 0.40 such that the ^-optimal solution 
is still guaranteed and also the converged speed is reasonably fast. 
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TABLE V 

Selection of e on the Performance of the Proposed Algorithm 



Value of e 


Number of iterations 


Converged WSR 


0.05 


8183 


4.7625 


0.10 


3498 


4.7438 


0.15 


2212 


4.7275 


0.20 


1642 


4.6942 


0.25 


1396 


4.6825 


0.30 


1148 


4.6620 


0.35 


1029 


4.6350 


0.40 


866 


4.6165 


0.45 


651 


4.5880 



SIMO Interference Channel 



- — Algorithm I 

* — Price-Based Algorithm 



20 40 60 80 100 120 140 
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Fig. 5. Performance comparison for Algorithm U versus the price-based algorithm in SIMO-IC. 

B. Providing Performance Benchmark for Other Heuristic Algorithms 

A key application of our proposed algorithm is to provide performance benchmarks for other heuristic algorithms 
for achieving the maximum WSR in the GIC, especially in cases of MISO-IC and SIMO-IC where the globally 
optimal solution by exhaustive search is hardly possible. In the following, we provide an example to show how to 
utilize our proposed algorithm to evaluate the performance of other suboptimal algorithms for WSR maximization 
in MISO-IC and SIMO-IC. 

We consider the "price-based" suboptimal algorithm, e.g., the ADP algorithm, which was proposed in (51 as 
an efficient distributed algorithm for WSR maximization in SISO-IC. Since to our best knowledge, extensions of 
ADP to the multi-antenna GIC are not yet available in the literature, we provide the details for such extensions 
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Fig. 6. Performance comparison for Algorithm U versus the price-based algorithm in MISO-IC. 



for SIMO-IC and MISO-IC in Appendix E 

Figs. [5] and [6] show the achievable sum-rates by the price-based algorithm versus Algorithm U for 4-user SIMO-IC 
and MISO-IC, respectively, without the minimum rate constraints. Each element in all channel vectors involved 
is randomly generated by the CSCG distribution with zero mean and unit variance. We set the parameters to 
control the accuracy of Algorithm U as e = 0.01 and 77 = 0.5. In Fig. \5\ the price-based algorithm converges 
to the sum-rate of 10.6989 in SIMO-IC, while the maximum sum-rate achieved by Algorithm U is 11.9182. In 
Fig. [6l the price-based algorithm converges to the sum-rate of 4.8216 (although it has reached almost 6 before 
convergence) in MISO-IC, while Algorithm U achieves the maximum sum-rate of 10.6193. Based on these results as 
well as other numerical examples (not shown in this paper due to the space limitation), we infer that in general the 
price-based algorithm for SIMO-IC performs better than MISO-IC, as both compared with our proposed algorithm 
that achieves the global sum-rate maximum. Moreover, the price-based algorithm for MISO-IC does not converge 
under certain channel setups, while even when the algorithm converges, the resulted sum-rate can be far from the 
global maximum. In contrast, for SIMO-IC, the price-based algorithm usually achieves the sum-rate very close to 
the global maximum, and even converges to it under certain channel setups. 

VI. Concluding Remarks 

In this paper, we propose a new framework to achieve the global optimality of WSR maximization problems in 
SISO-IC, SIMO-IC, and MISO-IC, with the interference treated as Gaussian noise. Although the studied problems 
are non-convex with respect to the power allocation and/or beamforming vectors, we show that they belong to the 
monotonic optimization over a normal set by reformulating them as maximizing the WSR in the achievable rate 
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regions directly. Therefore, the outer polyblock approximation algorithm can be applied to achieve the global WSR 
maximum. Furthermore, by utilizing the approach of rate profile, at each iteration of the proposed algorithm, the 
updated intersection point on the Pareto boundary of the achievable rate region is efficiently obtained via solving 
a sequence of SINR feasibility problems. It is worth noting that although the developed framework in this paper 
is aimed to solve the WSR maximization problem for the GIC, it can be similarly applied to other multiuser 
communication systems with non-convex rate regions provided that the problem of characterizing the intersection 
Pareto boundary point with an arbitrary rate -profile vector can be efficiently solved. 

It is worth pointing out that based on our numerical experiments, the proposed algorithm in this paper is found 
to converge very slowly when the number of users becomes large (e.g., K > 6), and thus may not be suitable for 
real-time implementation. Nevertheless, the proposed algorithm can be applied to provide performance benchmarks 
for other real-time algorithms that usually guarantee only suboptimal solutions. It is our hope that this paper will 
motivate future work to improve the convergence speed of the proposed algorithm and thus make it more applicable 
in practical systems, even with large number of users. For example, in 11331 , the original point for the algorithm 
is shifted from the origin to a point in the negative plane, which is shown to speed up the convergence to some 
extent. 

After the submission of this manuscript, we become aware of one interesting related work ||34l that is worth 
mentioning. In 11341 , a similar framework is proposed to optimize the system performance for multi-cell downlink 
MISO beamforming (similar to MISO-IC in nature), e.g., sum-rate performance and proportional fairness, by 
making use of monotonic optimization and rate profile techniques. One difference between [34] and our work 
is that for the monotonic optimization part, a so-called "branch-reduce-and-bound" algorithm is used in P4l as 
compared to the outer polyblock approximation algorithm in our paper. 

Appendix 

A. Proof of Lemma I4.il 

Note that under the sum power constraint, a similar result to this lemma has been shown in ||28l . However, the 
proof in ll28l is not directly applicable in our case since in (l36l ). there is an individual power constraint rather than 
the sum power constraint. Thus, we need to provide a new proof for this lemma shown as follows. 

Suppose that there are two solutions to (|36l ), denoted by (p,C(W)) and (p' ,C'(W)). Define a sequence of 
#fc's as Ok = j=JS V/c. We can re-arrange 9k 's in a decreasing order by 

e n >e t2 >--->d tK . (52) 

Since according to (l32l we have pi = p ■ = P™, it follows that 9i = 1 must hold. Hence, 9 t \ > 9% = 1. Moreover, 
in d52l , at least one strict inequality must hold because otherwise 9k = 1, Vfc, which then implies that only one 
unique solution to d36l ) exists. 
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Next, we derive the SINR balancing value of user tl as follows: 

P , tl \\w%hti,ti\\ 2 



C' tl (W) 



w a( E Pj h ti,jhn,j + <JtiI)wtiiti 

mi 

ptiWw^ht^tiW 2 

«>ti( E Pjhtijhgjfc + all-±)w n jti 

m||wg/i t i if i|| 2 

E Pjhtijh&j + a^wajti 
= C tl (W). (53) 

Based on (|3TT) . we have 

C'(W") = C^W) >C n {W) = C(W). (54) 
Similarly, we can show that C' tK (W) < CtK(W), which yields 

C'(W) = C' tK {W) < C tK (W) = C(W). (55) 
Since d54l and d55l ) contradict to each other, there must be one unique solution to ( f36b . Lemma |4~T1 is thus proven. 

B. Price-Based Algorithm for SIMO-IC and MISO-IC 

In this part, we provide the details of the suboptimal price-based algorithms for Problems (PI. 2) in SIMO-IC 
and (PI. 3) in MISO-IC, which can be viewed as extensions of the ADP algorithm proposed in [5] for SISO-IC. 
In ADP, each user announces a price that reflects its sensitivity to the interference from all other users, and then 
updates its transmit power by maximizing its own utility offset by the sum interference price received from all the 
other users. It was shown in Q that ADP can converge to the solution that has the same Karush-Kuhn-Tucker 
(KKT) conditions as that of the WSR maximization problem, and is thus guaranteed to achieve at least a locally 
optimal solution. In the following, we extend the ADP algorithm in [5] to SIMO-IC and MISO-IC, but without 
the proof of convergence. 

1 ) Price-Based Algorithm for SIMO-IC: 

In this part, we extend the ADP or price-based algorithm to SIMO-IC. First, without loss of generality, we 
substitute the optimal MMSE-based receive beamforming vectors for w^s into the SINR expression given in (0]). 
Then, given any transmit power vector p, the achievable rate for user k can be expressed as 

Rkip) = log 2 (l + 7fc SIMO - IC ) =log 2 (l+puh^Y^Pjh^h^ + alir^X (56) 

Thus in Problem (PI. 2), we only need to find the optimal transmit power solution, without the need to consider 
the receive beamforming optimization. 

Next, we present the KKT optimality conditions of Problem (PI. 2) with the objective function specified in (l56l ). 
For any locally optimal power solution p* , there exist unique Lagrangian multipliers A = (Ai, • • • , Xk) such that 



25 



for any k = 1, • • ■ , K, 

dR k (p*) ^ dRjip*) 



Pk- 



dpk 

3+k 



X k (Pr x -Pk) = 0, (58) 
A fc > 0. (59) 

Now, for the price-based algorithm, define the price charged by receiver j to transmitter k, which indicates the 
sensitivity of the achievable rate of receiver j subject to the power change of transmitter k, as 

aD , , p., i'" ay. Pihi,ih?,i + °]ir%,k\\ 2 
- = — i v (60) 

In 2 ( 1 + pjhf^Epihj,^ + o)l)- l h hi \ 

Consequently, we see that the KKT conditions in (1571 ), (|58l ) and d59l ) are both necessary and sufficient for the 
optimal solution to the following problem for user k, k = 1, • • • , K : 

Maximize p k log 2 ( 1 + p k h kk (^2 i p j h k jh kj + oll)~ l h kM j - p k ^ Pj^j,k 

Subject to p fc < P fc max , (61) 

where pj and vr^fc are fixed, \/j ^ k. 

Similar to the ADP algorithm in 0, we propose the following algorithm to update the price and transmit power 
iteratively for all users in SIMO-IC. Specifically, at each iteration the algorithm does the following: 

1. Each user announces its price obtained using (f60b to all the other users; 

2. Each user updates its transmit power by solving Problem (l6Tb . i.e., 



Pk 



/'/,• 1 



In 2 £ m , k h? k { £ Pjhkjh^j + a 2 k I)- l h Kk 



, VA;, (62) 
o 



where [x]^ = max(min(x, b), a). 

Because Problems (PI. 2) and (f6Tb possess the same KKT optimality conditions, when the above algorithm 
converges to a set of optimal solutions to problems in (f6Tb for all fc's, this set of solutions will be at least a locally 
optimal solution to Problem (PI. 2). 

2 ) Price-Based Algorithm for MISO-IC: 

Next, we extend the ADP algorithm to MISO-IC. For any given transmit beamforming vectors V, we first define 
the price for user k as 

= _dRk = ll^fcll 2 

* k dT k ln2(\\h» k v k \\2+r k + al)(r k + aiy 

where F k = ^ \\h k jVj\\ 2 is the total interference power at the kth receiver. Let S k = v k v k , \/k. Given fixed 

interference prices and beamforming vectors of all other users, the following problem is to be solved by any user 
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k for its own transmit beamforming update: 

Maximize fi k log 2 (l + k ' k * |^ ) - ^ HjKjh ktj S k h jik 

Subject to Tr(5 fe )<P fc max 

S k y 0, (64) 

where S k >z means that is a positive semi-definite matrix. Similar to the previous case of SIMO-IC, we can 
show that the KKT conditions of Problem (l64l with k = 1, • • ■ , K are also those of Problem (PI. 3) by replacing 
VkV k with Sk, V/c. However, Problem (PI. 3) requires that the optimal solution S* k in Problem (l64l to be rank-one, 
which is not guaranteed a priori. Thus, Problem d64l is a relaxation of the original WSR maximization problem 
(PI. 3) without considering the rank-one constraint. 

Interestingly, it was recently shown in ifTOll that the optimal solution to Problem (l64l is always of rank-one, 
i.e., S k = v k v k . Hence, we propose a price-based algorithm for MISO-IC in a similar way to that for SIMO-IC. 
When this algorithm converges to a set of optimal solutions to problems in (l64l with k = 1, ••• ,K, this set of 
solutions are all rank-one and thus corresponds to at least a locally optimal solution to Problem (PI. 3). 

For this price-based algorithm for MISO-IC, the interference price can be iteratively updated according to (l63l . 
As for the update of beamforming vectors, we need to solve Problem (l64l for each user k. It can be verified that 
Problem (l64l is convex with strictly feasible points, and thus it can be solved by the standard Lagrangian duality 
method ll24l with a zero duality gap. The details of solving Problem (l64b can be found in Appendix I of iTTOl . and 
are thus omitted here. 
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